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Abstract 

In this paper, we address the problem of globally stabilizing a linear time-invariant 
(LTI) system by means of a statie feedbaek law whose amplitude and sueeessive 
time derivatives, up to a preseribed order p, are bounded by arbitrary preseribed 
values. We solve this problem for two elasses of LTI systems, namely integrator 
ehains and skew-symmetrie systems with single input. For the integrator ehains, 
the solution we propose is based on the nested saturations introdueed by A.R. Teel. 
We show that this eonstruetion fails for skew-symmetrie systems and propose an 
alternative feedbaek law. We illustrate these findings by the stabilization of the 
third order integrator with preseribed bounds on the feedbaek and its first two 
derivatives, and similarly for the harmonie oseillator with preseribed bounds on 
the feedbaek and its first derivative. 

Keywords: Global stabilization. Bounded eontrols. Rate saturation. Multiple 
integrators, Oseillators. 


1. Introduction 

Actuator constraints constitute an important practical issue in control appli¬ 
cations since they are a possible source of instability or performance degrada- 
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tion. Strong research efforts have been devoted to the stabilization of linear time- 
invariant (LTI) plants. LTI systems are known to be global stabilizable despite 
actuator saturations (i.e., by bounded inputs) if and only if they are stabilizable 
in the absence of input constraints and their internal dynamics has no eigenvalues 
with positive real part [[ll]. For systems that do not fulfill these constraints, several 
approaches provide stabilization from an arbitrarily large compact set of initial 
conditions (semiglobal stability); these include for instance ylsIdtlSD. Some of 
these semiglobal approaches can be extended to robust stabilization in presence 
of exogenous disturbances 0,0]. 

The objective of globally stabilizing plants by a bounded feedback remains of 
practical relevance, since the resulting control gains do not depend on the magni¬ 
tude of initial states. Procedures ensuring global stability of nonlinear plants by 
bounded feedback have been proposed in 0 0] and robustness to exogenous in¬ 
puts have been addressed in liloL llin. Among the LTI systems that can be globally 
stabilized by bounded feedback, chains of integrators have received specific atten¬ 
tion. The simple saturation of a linear feedback fails at ensuring global stability as 
soon as the integrator chain is of dimension greater than or equal to three [10, 131. 
In llldll a globally stabilizing feedback was constructed using nested saturations 
for a chain of integrators of arbitrary length. This construction has been extended 
to all LTI plants that can be stabilized by bounded feedback in [HI], in which a fam¬ 
ily of stabilizing feedback laws was proposed as a linear combination of saturation 
functions. In (IlSll . the issue of performance of these bounded feedbacks is inves¬ 
tigated for chains of integrators and some improvements are achieved by using 
variable levels of saturation. Global approaches ensuring robustness to exoge¬ 
nous disturbances have also been investigated. The first general solution to the Lp 
finite-gain stabilization problem was provided in based on a gain scheduled 
feedback initially proposed in lll7n . An alternative easily implementable solution 
to that problem was recently proposed in [[101 for chains of integrators. As for 
neutrally stable systems, such a solution was first given in s. 

While actuation magnitude is often the main concern in practical applications, 
limited actuation reactivity can also be an issue. Indeed, technological constraints 
may affect not only the amplitude of the delivered control signal, but also the 
amplitude of its time derivative. This latter problem is known as rate saturation 
and has been addressed for instance in EOl l2ll . 1221 l23l . 12411 Semiglobal stabiliza¬ 
tion has been achieved via a gain scheduling techn iq ue EQII . or through low-gain 
feedback or low-and-high-gain feedback [24]. In [20 23ll. regional stability was 
ensured through LMI-based conditions. In 112ill , this problem has been addressed 
for nonlinear plants using backstepping procedure ensuring global stability. 
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In this paper, we deepen the investigations on global stabilization of integra¬ 
tor ehains subjeet to bounded aetuation with rate eonstraints. We eonsider rate 
eonstraints that affeet not only the first time derivative of the eontrol signal, but 
also its sueeessive p first derivatives, where p denotes an arbitrary positive inte¬ 
ger. We speeifieally study two elasses of systems that ean be globally stabilized 
by bounded state feedbaek, namely ehains of integrator and skew-symmetrie dy- 
namies with single input. No restrietion is imposed on the dimension of these 
systems. We show that solving the problem for these two oases aotually oover 
wider elasses of systems, namely all systems with either only zero eigenvalues or 
only simple eigenvalues with zero real part. For both these elasses of LTI systems, 
we propose a bounded statio feedbaek law that ensures global asymptotio stability 
of the olosed-loop system, and whose magnitude and p first time derivatives are 
bounded by arbitrary preseribed values. For the ehains of integrators, the pro¬ 
posed eontrol law is based on the nested saturations prooedure introduoed in [Il4l] . 
We rely on speeifie saturation funetions, whieh are linear in a neighborhood of the 
origin and eonstant for large values of their argument. Unfortunately, we show 
that this nested saturations feedbaek fails solve the problem for skew-symmetrie 
dynamies. For the latter elass of systems, we propose an alternative eonstruetion. 

This paper is organized as follows. In Seetion [2l we provide definitions and 
state our main results for both eonsidered elasses of LTI systems. The proofs 
of the main results are provided in Seetion [phased on several teehnieal lemmas. 
In Seetion S we test the effieieney of the proposed eontrol laws via numerieal 
simulations on the third order integrator and the harmonie oseillator, where we 
bound with preseribed values the feedbaek, as well as its first two time derivatives 
for the third order integrator and it first time derivative for the harmonie oseillator 
respeetively. 

Notations. The funetion sign : M\{0} —)■ R is defined as sign(r) := r/|r|. 
Given a set 7 C R and a eonstant a G R, we let I>a '■= {vG / : x> a}. Given 
k G N and n,p G N>i, we say that a funetion / : R" —)■ R^ is of elass C*(R",R^) 
if its differentials up to order k exist and are eontinuous, and we use to denote 
the k-th order differential of /. By eonvention, := /. The faetorial of k is 
denoted by k\ and the binomial eoeffieient is denoted (^) := —rrP —rr- We define 
|m,k] := {/ G N : / G [m,k]}. We use R"’” to denote the set of n x n matriees 
with real eoeffieients. The matriees In and G R”’" denote the identity matrix of 
dimension n and the n-th Jordan bloek respeetively, i.e., the nxn matrix given by 
{Jn)ij = 1 if i = 7 — 1 and zero otherwise. For eaeh i G |1, n], e,- G R" refers to the 
eolumn veetor with eoordinates equal to zero exeept the i-th one equal to one. We 
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use ||a:|| to denote the Euelidean norm of an arbitrary veetor x G M”. For two sets 
A and B, the relative eomplement of A in 5 is denoted by B\A. 


2. Statement of the main results 

2.1. Problem statement 

We start by introdueing in more details the general problem we address. Given 
n G N>i, eonsider LTI systems with single input: 

x = Ax-\-Bu^ (1) 


where a: G M", A and B are nxn and n x 1 matriees respeetively. Assume that 
the pair (A,5) is stabilizable and that all the eigenvalues of A have non positive 
real parts. Reeall that these assumptions on (A,5) are neeessary and suffieient 
for the existenee of a bounded eontinuous state feedbaek u = k{x) whieh globally 
asymptotieally stabilizes the elosed-loop system [U]. We say that an eigenvalue of 
A is critical if it has zero real part. 

Given a family of preseribed bounds {Rj)o<j<p on the eontrol signal and its 
sueeessive p-first derivative, we start by introdueing the notion of p-bounded feed¬ 
baek law by {Rj)o<j<p for system ([T]). This terminology will be used all along the 
doeument. 


Definition 1. Given n G N>i and p G N, let {Rj)o<j<p denote a family of positive 
constants. We say that v : R" ^ M is a p-bounded feedback law by {Ri)o<j<pfor 
system ([7]) if it is of class and, for every trajectory of the closed-loop 

system x = Ax-\-Bv{x), the control signal U : M>o —)■ M defined by U{t) := v{x{t)) 
for allt>0 satisfies, for all y G |l,p], sup{|t/(-^^(t) | : t > 0} < Rj. 

Based on this definition, we ean restate our stabilization problem as follows. 


Problem. Given p G N and a family of positive real numbers {Rj)o<j<p, design 
a feedback law V such that the origin of the closed-loop system x = Ax^-By{x) is 
globally asymptotically stable (GAS for short) and the feedback v : M" ^ M is a 
p-bounded feedback law by {Rj)o<j<pfor System ([7]). 


The ease p = 0 eorresponds to ^bal stabilization with bounded state feedbaek 
and has been addressed in e.g. IJA.Ul]. The ease p = 1 eorresponds to global 
stabilization with bounded state feedbaek and limited rate, in the line of e.g. 


23Ll2d,lMl2ll]. 
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In this paper, we present a general solution to the problem at stake for two 
classes of LTI systems: Case 1 all the critical eigenvalues of A are zero, Case 2 
all the critical eigenvalues of A are simple and have zero real parts. 

Since the pair (A, 5) is stabilizable, there exists a linear change of coordinates 

transforming the matrices A and B into j and f ’ where A i is Hurwitz, 


the eigenvalues of A2 have zero real parts and the pair (A2,52) is controllable. 
Then, it is immediate to see that we only have to treat the case where A has only 
critical eigenvalues. From now on, we therefore assume that A has only eigenval¬ 
ues with zero real part, and that the pair (A, 5) is controllable. 


2.2. Multiple integrators 


In Case 1, up to a linear change of coordinates, A can be put in a block- 
diagonal form with Jordan blocks Jr on the diagonal. It is then clear that, up to 
an additional linear change of coordinates, addressing Case 1 amounts to dealing 
with the sole case of a multiple integrator of arbitrary length n, i.e. the LTI control 
system given by 


Letting X := (xi, 


i<n and x„ = u. (2) 

system Q can be compactly written as i = J„x + enU. 
Inspired by Hldli . our design of a p-bounded feedback for this system is based 
on a nested saturations feedback. We focus on the specific class of saturations that 
are linear around zero, and constant for large values of their argument. 


Definition 2. Given p G N, ^{p) is defined as the set of all odd functions o of 
class C^(M, M) such that there exist positive constants a, L, and S satisfying, 
for all r G M, ro{r) > 0 when r 7^ 0, o{r) = ar for all \r\ < L, and |o'(r)| = 
^max^ wfien |r| >S. In the sequel, we associate with every o E ^{p) the 4-tuple 



Figure 1: A typical example of a =5^(p) saturation function with constants ,L,S, a). 

The constants a will be widely used throughout the paper, see Fig. 

[Uto fix ideas. Notice that it necessarily holds that S>L and the equality may only 
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hold when p = 0. We also stress that the sueeessive derivatives up to order p of an 
element of ^ (p) are bounded. An example of sueh funetion is given in Seetion 
l4.2l for p = 2 . The first result of this paper establishes that global stabilization on 
any ehain of integrators by bounded feedbaek with eonstrained p first derivatives 
ean always be aehieved by a partieular ehoiee of nested saturations. In other 
words, it solves the Problem in Case 1. 

Theorem 1 . Given n G N>i and p E N, let {Rj)o<j<p be a family of positive 
constants. For every set of saturation functions 0'i,...,a„ G y{p), there exist 
vectors ,..., in M", and positive constants ai,..., such that the feedback 
law V defined, for each x G M", as 

v(v) = -anOn{klx + an-\On~\{kl_iX+ ... + a\0\{klx)) ... j ( 3 ) 

is a p-bounded feedback law by {Rj)o<j<p far system (|2l), and the origin of the 
closed-loop system x = JnX + e„v(x) is GAS. 

The proof of this result is given in Seetion ITT] and the argument also provides 
an explieit ehoiee of the gain veetors k\,.. .,kn and eonstants ai,... a„. 

Remark 1 . Along the proof of Theorem\T\that the proposed construction allows 
to chose the magnitude of control signal independently of the magnitude of its p 
first derivatives. More precisely, can be chosen to ensure that max{|v(x)| : x G 
M”} = Rq and the gain vectors ki,.. .,kn and constants ai,... a„_i can be taken in 
such a way that the p first derivatives of the feedback are bounded by (Rj) i<j<p- 

Remark 2 . In /Q/, a stabilizing feedback law was constructed using linear com¬ 
binations of saturated functions. That feedback with saturation functions in 5 A [p) 
cannot be a p-bounded feedback for System (|2])/or p>\. To see this, consider the 
double integrator, given by xi =X2,X2 = u. Any stabilizing feedback using a linear 
combination of saturation functions in 5 A (p) is given by v(xi ,^2) = —aO\ {bx2) — 
C02{d{x2 +xi)), where the constants a, b, c, and d are chosen to insure stability of 
the closed-loop system according to ft Let U (t) = v(xi (t),X2{t)) for all t > 0 . 
A straightforward computation yields, for t > 0, U{t) = —abo^\ax 2 {t))U[t) — 
cd02\d{x2{t) -\- x\{t))){x2{t) -\-U {t)). Consider now consider a solution with ini¬ 
tial condition ^2(0) = X20, and xi(0) = —X20 such that o[^\ax2o) = 0. Wb then 

have U{ 0 ) — —c<ia2^^(0)(x2o + C( 0 )), whose norm is greater than ci(|x2o| —C2) 
for some positive constants c\,C2. Thus |ii( 0 ) | grows unbounded as |v2o| tends to 
infinity, which contradicts the definition of a 1 -bounded feedback. 
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2.3. Harmonic oscillators 

In Case 2, up to a linear ehange of eoordinates, A can be put in a block- 
diagonal form with skew-symmetric matrices on the diagonal. Addressing the 
stabilization problem under concern therefore amounts to only considering the 
following control system : x = Ax-\~bu, where x G M”, A G M"’" is skew-symmetric, 
Z? G M” and the pair (A, b) is controllable. 

Unfortunately, the nested saturations feedback law given in ([3]) is a generic 
solution to the Problem for this class of systems only in the scalar case (n = 1) 
or when when no rate constraint is imposed (p = 0). To see why it may fail for 
n>2, consider for instance the harmonic oscillator given by ii =X 2 ,X 2 = —xi -t- u 
(which we address in more details in Section |4^ with a bounded stabilizing law 
given by M = —o{x2) with o G y{p) for some integer p. The time derivative 
of M then satisfies \u{t) \ > |a*^^)(v 2 (t))|(|.ri(t)| — |M(t)|), which grows unbounded 
as xi goes unbounded and 1^21 remains small (i.e. in the linear zone of a). This 
prevents the feedback —o{x 2 ) to be a 1-bounded feedback, hence a p-bounded 
feedback for all p > 1. 

Our second result provides an alternative p-bounded feedback for the har¬ 
monic oscillator, thus solving the Problem in Case 2. 

Theorem 2. Given n G N>i and p G N, let iRj)o<j<p be a family of positive 
constants, let A G M”’” be a skew-symmetric matrix, and let b ^W' be such that 
the pair {A,b) is controllable. Then, for any a > 1/2, there exists a positive 
constant /3 such that the feedback law v : M" —?■ M defined as 

v(v):=-/3hV(l + IWl")“) (4) 

is a p-bounded feedback law by {Rj)o<j<pfor x = Ax-\-bu, and the origin of the 
closed-loop system X =Ax-\-b'v{x) is GAS. 

The proof of this theorem is given in Section [X2l 

Remark 3. Unlike for multiple integrators (see Remark\T^, the magnitude of the 
proposed feedback is not independently of the amplitude of its p first derivatives. 

3. Proof of main results 

3.1. Multiple integrators 

We start by estimating upper bounds on composed saturation functions of the 
class ^(p). This estimate, presented in Lemma [2l relies on Faa di Bruno’s for¬ 
mula recalled below. 
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Lemma 1 (Faa Di Bruno’s formula, [l25n . p. 96). Given k eN, let (j) e C^(M>o,ffi) 
and p G C*(M, M). Then the k-th order derivative of the composite function p o ^ 
is given by 


|po<,lW(,)= f (5) 

a=\ 

where the Bell polynomial given by 

8e^k,a l=l 

where denotes the set of{k — a+ \ )—tuples 8 := (5i, § 2 , ■ ■ •, ^k-a+i) of pos¬ 
itive integers satisfying 5i + 62 + ... + 8k-a+i = o and 5i + 2^2 + {k — a-{- 

l)4-fl+i =k, andc§ := k\/ (^5i\■ ■ ■ Sk-a+iK'^'Y^ ■ ■ ■ {{k-a+l)\)^’<-a+iy 

Remark 4. We stress that the Bell polynomial B^ a of (homogeneous) degree a 
w.r.t the (k — a-\- \ )-dimensional vector representing the argument ofB^ a- 

The proof of Theorem [T] extensively relies on the following upper bound on 
composition of functions of S^{p), which exploits their constant value in their 
saturation region. 

Lemma 2. Given G N, let f and g be functions of class C*(M>o,M), o be a 
saturation function in y (k) with constants ( L, S, a), and E and F be subsets 
o/M>o such that EOF. Assume that 

|/(t)|>5, VtGn^, (7) 

and that there exist positive constants M,Q\,... ,Qk such that 

|g(^)(OI<M,VtGF and |/(“)(0| < <2a, Vt G Va G |l,fc]. (8) 

Set 'aa '■= max{|o'(‘^)(5')| : 5 G M} for each a G |l,fc]. Then the kth-order 
derivative of the function h : M>o M defined by h(-) := g(-) + o(f(-)), satis¬ 
fies 

k 

\hf)(t) \<M+Y, aaBk,a{Qu- ■ ■, Qk-a+l): Vt G F. (9) 

a=\ 
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Proof of Lemma^ Using Lemma [B a straightforward computation yields 

*“'(») = v<> 0 , 

a=\ 

where the polynomials 5^^ are defined in ®. Since c G © ensures that 

the set F \£' is contained in the saturation zone of a. It follows that 

[cyo/]W(t) =0, VteF\E. (10) 

Furthermore, from ® and ® it holds that, for all t G F, 


k—a+l 


C5 n Q?=Bk,a{Ql^---^Qk-a+l)- 
/=1 


( 11 ) 

One has [c 7 o/](^)(t) <Y!'a=l(^aBk,aiQl^---^Qk-a+l), for all t G F, from the defi¬ 
nition of and ®. In view of (fTO . the above estimate is valid on the whole set 
F. Thanks to ®, a straightforward computation leads to the estimate dH). □ 


We now turn to the proof of Theorem [B We explicitly construct the vec¬ 
tors and the constants guaranteeing global asymptotic sta¬ 

bility with a bounded feedback law whose successive derivatives remain below 
prescribed bounds at all times. Given p G N and n G N>i, let a,- be saturation 
functions in y'{p) with constants ttcj,.) for each i G |l,n], and let 

(Fy)o<y<p be the family of prescribed positive bounds on the amplitude and the 
successive time derivatives of the control signal. We first construct a p bounded 
feedback by {Rj)o<j<p for System dB). Then we show that the origin of the closed 
loop system dB) with this feedback law is GAS. 

Let (/i,™^)i<Kn-i and (F^,.)i<,<„_i be two families of positive constants 
such that < l/2,L^„_j = o'lff and, for each i e |l,n-2], 

pfax <; j2,Lp. = Ufj ./. For each fG |l,n — l],we define satura- 

tion function m eY(p) with constants (u™^, L„., Su,, 1), where Sm. = Sa,L„./La:, 
as follows pfs) := At™“cJ/(5La,./F^.)/a'"“, for all 5 G M. For A > 1, to be cho¬ 
sen later, we define the saturation function p„ as Pn{s) ■ = Ro(7„{sL(y,JX,) 
for all 5 G M, with jU/f" = Rq, = A, = Sa„^/La„, and := ag/A with 
ag :=FoLa„acT„/tyr^- 

We next make a linear change of coordinates y = Hx, with H G M"’", that puts 
System dB) into the form 


3),- = ag/A yi + u, V/G|l,n]. 

l=i+\ 


( 12 ) 
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The relations yn-i = T.k=o {[) for i ^ [0,^ “ ll> enable us to deter¬ 

mine H. Since H is triangular with non zero elements on the main diagonal, it is 
invertible. We define a nested saturations feedback law T : M” —)■ M as 

Tf(3^) = -l^n{yn + l^n-l{yn-l+--- + l^l{yi))---)- (13) 

Note that, in the original v-coordinates, this feedback law reads T(>') = T{Hx) 
therefore the bounds of the successive time derivatives of T(>') coincide with those 
of T (Hx). The global stabilization of (fT^ with a /^-bounded feedback law by 
{Rj)o<j<p is thus equivalent to that of the original system So, from now on, 
we will rely on this expression. Let y(-) be a trajectory of the system 

yi = ap/X ^ y/+T(y), Vie|l,n], (14) 

l=i+\ 

which is the closed-loop system (O with the feedback defined in (fTSl) . For each 
i e |l,n], let Zi : M>o —)■ M be the time function defined recursively as Zi{-) := 
y,(-) -t- /i,_i(z,_i(-)) , with /d(-) = 0. With the above functions, the closed-loop 
system (fT4l) can be rewritten as 

{ n—i 

y/ = ag/Az„-Atn(zn) + ag/A( E yi{zi)-^i{zi)),\/iell,n-l\, 

l=i+\ (15) 

yn — l^n{Zn)i 

where fiizi) — Zi — IJ-i{zi) for all I G |l,n — 1]. For each i G |l,n], we also let 
Ei := {y G M" : |yv| < +G and /,■ := {t G M>o : y(0 G 

Note that from the definitions of 7, and we have /i C 72 C ... C 7„, and a 
straightforward computation yields 

|z,'(t)| Vt G 7/+i\7/, Vz G |l,n-1], (16) 

\Zn{t)\>Sp„, VtGM>o\7„, (17) 

which allows to determine when saturation occurs. We define bp. := max{|r — 

l^i{f)\ '■ \f\ + and/I,-^^-:=max{|/ip^(r)|: r G M} for each z G |l,zz —1] 

and j e |l,p]. Note that these quantities are well defined since the functions /i,- 
are all in y'{p). We can now establish that 

\zi{t)-iJ,i{zi{t))\<bp., Wteli, VzG|1,zz-1]. (18) 
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Set := min{(7„(r)/r : 0 < |r| < Sf 7 „ +2jU™fLa„}, ba„ := max{o-,j(r)/r : Vr > 
0}, and A = {bo,, —bfjJ{Lo„Ro/(yn‘‘^)- It then ean been seen that 

\a^,„Zn{t)-^n{Zn{t))\ < A(Sct„/L,,„ + 2^^A), Vt G 4- (19) 

The following statement provides explieit bounds on the sueeessive derivatives 
of each functions yi{t), Zi{t) for each i G |l,n] and the control input given by 

['(■):= TWO). 

Assertion 1. With the notation introduced previously and the Bell polynomials 
introduced in ®, every trajectory of the closed-loop system f[7?l) satisfies, for 
each i G |l,n] and each j G |l,p], 

Wi(i'J)): bF’WIST.j.Vte/i; {PiiiJ))- kF’WISZij.V'e/i, (20) 
(P,(j)) : sup{|(/((l(»)| :! > 0} < ^ (21) 

q=i 


where = max{|/ri^^(r)| : r G M}, Yij, Zij, and Gqj are independent of initial 
conditions and are obtained recursively as follows: for j = 1, T) i := A{So„/Lo„ + 
2p'ffj/X)-Ya^/HL’;zU^U,+l^r)forielfn-llYn,, :=pr^ Zi,i -Tpi, 
:= y,-,i for i G | 2 ,n], and Gyi := Z„j. For each j G | 2 ,p], 


the inductive relations are given by Yij := (Xfi/^Y!b=i+i ^bj~i +L^=i 

for i E |l,n], Zyy := Yij, Zjj := Yij + Y,a=i j-_i 4 _a),/or 

iE |2,n], andGqj:=Bj^q{Zny,...,Znj-q+i),forqE |1,7]. 


Proof of Assertion\T\ The right-hand side of System (fTO is globally Lipschitz and 
of class CP{'MA,'MA), and therefore it is forward complete with trajectories of class 
C^+^(R>o,M"). We establish the result by induction on j. We start by j = 1. 
We begin to prove that Pi{i, 1) holds for all i E |l,n]. Let i E |l,n — 1]. From 
(fTSl) . (fT^ . and (fT9l) a straightforward computation leads to |yi(t)| < A{S(y„/L(y„ -+- 
-f for all t E //+i. Since f C the above 

estimate is still true on /,■. Moreover, from (fTSl) it holds that \yn{t) \ < at all 
positive time. Thus, Pi{i, 1) has been proven for each i E |l,n]. 

We now prove by induction on i the statement P 2 {i, 1). Since zi (■) = yi (•)> the 
case i = 1 is done. Assume that, for i E |l,n — 1], Piib 1) holds. From Lemma 
|2](with 1, /(•) ^Zi{-), g{-) =yi+i{-), hf) =z,+i(.), a(-) = pf-), Qi =Zi^i, 
M = 1 , ai =/!■ J, £ = Ii, F = 4+1, and (fT^l. we canestablish thatF 2 (i +1,1) 

holds. Thus P 2 {i, 1) holds for all z G |l,zz]. 
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Notice that the applied eontrol input reads U{-) = —jUn(zn(-))- We then ean 
establish ^3(1) from Lemma [2] (with k = 1, /(•) = z«(-), g(-) = 0, /?(•) = f/(-), 
o-(-) = Qi = Z„j, M = Q,Oi= ;tr„ 1 , E = In, F = R>o and (QT])). This ends 
the ease 7 = 1. 

Now, assume that for a given j e |l,p — 1 ], statements P\{i,j2), Piiiji) 
and P'i{i2) hold for all 72 < 7 and all i G |l,n]. Let i E |l,n]. From (fl^ . a 
straightforward eomputation yields < ocp,/^,^4=1+1 

for all t > 0 . Using ^3(7), Fi(/+ 1,7),... ,Fi(n,7), we obtain that |yp^^^( 0 l < 
CCji/X L'/=,+i T/,/ + Yjq=i Gq,jPn,q ’ f^r all It > li. Thus the statement Pi (7 + 1 , i) is 
proven for all z G |1 , zz]. 

We now prove by induetion on i the statement ^ 2(^7 + 1). As before, sinee 
zi(-) = yi(-)> the ease for z = 1 is done. Assume that for a given z G |l,zz — 1], 
the statement F 2 (z ,7 + 1) holds. From Lemma |2] (with k = 7 + 1, /(•) — Zi{-), 

^(•) = K-) = Z/+l(-)> C^(-) = Qh = ZiM’M = Yi+Ij+I, Oa = Pi^a, 

£ = 4F=4+1, and (fT^l. we ean establish that/’ 2 (z +1,7 + 1) holds. P 2 {i,i+l) 
is thus satisfied for all z G |1, zz]. 

Finally, we ean establish ^3 ( 7 + 1) from Lemma[2]( with = 7 + l,/(-) —Zn{-), 
g(-) = 0, /z(-) = [/(■), C7(-) = jUn(-), Qkl = Zn,ki, M = 0,aa= Jin,a^ E = !„, F = 
M >0 and (fT7]) l. This ends the proof of Assertionfl] □ 


First notiee that 7 t„ ^ = iln,qlZ^ with jln,q := and Gn,q := 

max{|o'(-^)(r)| : r G M} for ^ G |l,p|. In eonsequenee, using AssertionfH it ean 
be seen that, for eaeh 7 G |l,p|, Y.q=iGqjjl„^q/= F(l/A)/A where P is a 
polynomial with positive eoeffieients. This sum is thus deereasing in A. Henee, 
we ean piek A > 1 in sueh a way that Eg=i Gqjjln,q/X‘^ < R, for each 7 G |l,p| 

withP := minjPi,. ..,Rp}. It follows that, for each 7 G |l,p|, sup{|t/('^^(t)| : t > 
0} < P < Rj. By recalling that the feedbaek T is bounded by Rq, we eonelude that 
it is p-bounded feedbaek law by (P 7 )o< 7 <f- for System (fT^ . 

It remains to prove that the feedbaek law (fT3l) . where now all the eoeffieient 
have been ehosen, stabili z es System (fT^ . Aetually the proof is almost the same 
as the one given in [Il4ll . exeept that we allow the first level of saturation /i„ to 
have a slope different from 1 . 

We prove that after a finite time any trajeetory of the elosed-loop system (fT4l) 
enters a region in whieh the feedbaek (flSl) beeomes simply linear. To that end, we 
eonsider the Lyapunov funetion eandidate V„{y„) := Its derivative along the 
trajeetories of (fT4l) reads y„ = —ynPn {jn + Pn -1 {zn- 1 )). From the ehoice of A and 
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, we obtain the following implieation, with 6 : 


inf 




\yn\>Lfij2 %<-dL^j2. 


(22) 


We next show that there exists a time Ti > 0 such that \y„(t) \ < Lfi„l2, for all 
t > Ti. To prove that, we consider the following alternative: either for every t > 0, 

1 — ^Mn/2 we are done, or there exist Tq > 0 such that |y„(ro)| > L^J2. 
In that case there exists Tq > Tq such that yn(fo) = T^„/2 (otherwise thanks to 
(l22l) . Vn(t) as t which is impossible). Due to (l22l) . we have |>'n(t)| < 

L^,y2 in a right open neighbourhood of Tq. Suppose that there exists a positive 
time fi > Tq such that |y„(7'i)| > T^„/2. Then by continuity, there must exists 
fi e (Tq^Ti] such that \y„(f 2 ) \ = L^J2, and |y„(t)| < L^J2 for all t G (TqJ 2 ). 
However, it then follows from (l22l) that for a left open neighbourhood of T 2 we 
have \yn{t) \ > bn(? 2 ) I = T^,J2. This is a contradiction with the fact that on aright 
open neighbourhood of Tq we have |y,i(0l < Tjx,j2. Therefore, for every fi > Tq, 
one has |yn(ri) | < L^„/2 and the claim is proved. 

In consequence we have that |y„(t) + jU,i-i(z«-i(t))| < for all t >Ti. 
Therefore /i„ operates in its linear region after time Ti. Similarly, we now con¬ 
sider Vn -1 (y?!- 1 ): = \y'i- 1 , whose derivative along the trajectories of (O satisfies 
K = -1 Mn-1 (y «-1 + Mn -2 (yn -2 + ...)), for all t > Ti. Reasoning as before, 
there exists a time T 2 > 0 such that |y„_i(t)| < and /r„_i operates in its 

linear region for all t >T 2 . 

By repeating this procedure, we construct a time Tn such that for every t >Tn 
the whole feedback law becomes linear, that is T(y(t)) = —cc^„(yn(t) + ... + 
yi (t)), when t>T„. Thus, after time Tn, the system (fT4l) becomes linear and its lo¬ 
cal exponential stability follows readily. Thus the origin of the closed-loop system 
(fT4l) is globally asymptotically stable. 

With the linear change y = Hx proposed in (fT4l) . the closed-loop system (fT4l) 
can be put into the form of the original system Q in closed loop with u = 
T(Hx). Thus, the sought feedback law v of Theorem [T] is obtained by v(x) = 
T(Hx). This leads to the following choices of parameters: for each z G |l,n — 1, 



3.2. Skew-symmetric systems with scalar input 

Given n G N>i, let A G M”’” be a skew-symmetric matrix and b G M" such that 
(A,b) is controllable. Let p G N, a > 1/2 and (Rj)o</<p be a family of positive 
constants. The proof of Theorem |2] is divided into two steps. We first prove that 
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for any /3 > 0 and a > 1/2, the origin of x = Ax-\-bu with the feedbaek law © 
is GAS. We then show that for any a > 1/2 there exists a positive eonstant /3 
sueh the feedbaek law @ is a p-bounded feedbaek law by {Rj)Q<j<p for system 
X = Ax + bu. 

Let /3 be a positive eonstant and a > 1/2. We define Ap := A — f^bb^. The 
matrix A |3 is then Hurwitz. To see this, observe that the Lyapunov equation + 
Ap = —2^bb^ holds and that the pair (Ap,b) is eontrollable. Thus there exists a 
real symmetrie positive definite matrix P E M”’” sueh that PAp +A^P ~ —/„. 

Let V{x) := x^PxK{{ \ + ||x||^)“+^ — 1) be a eandidate Lyapunov funetion 
with K := {\\Pb\\ /{a + 1). The derivative of V along trajeetories ofx= Ax+ 

bv{x) is then given by 

V{x) = - ||xf+ 2/3(1 - 1/(1 + \\xf)^)x'^Pbb'^x-2K{a + l){b'^x)^. 

Using that \x'^Pb\ < \\x\\ \\Pb\\ and 

2/3(1-1/(1+ ||xf)“)|x^PZ7Z>^x| < \\xf/2 + 2[5^\\Pbf{b^x)^, 

we get that V(x) < — ||x|| /2. Therefore the origin of i = Ax + bv{x) is GAS. 

We next give and prove two assertions whieh give explieit formula of the sue- 
eessive derivatives of the trajeetories and of the eontrol signal. 

Assertion 2. Given any /3 > 0 and a > 1/2, each trajectory x : M>o —)■ K" of 
X = Ax-\-bv{x) is C°° and satisfies, for any k E Pl>i, with Uf) := v(x(-)), 


x(^)(0 =A’^x{t) + Vt > 0. 

;=o 


(23) 


Proof of Assertion^ The right hand side of i = Ax + bv{x) is globally and Lip- 
sehitz C°°, and therefore this system is forward eomplete with trajeetories are of 
elass C°° (M>o, ffi”). In partieular, the sueeessive derivatives of t/(-) and x(-) are 
well defined. Equation (l23l) then follows by a trivial induetion argument using 


differentiation of x = Ax + Z?v(x) at any order. 


□ 


Assertion 3. Given /3 > 0 and a > \/2, let x : M>o —)■ M” be any trajectory of 
x=Ax-\-bv{x) and let G(-) := 1 + ||x(-)||^. Then, for any k E N>i and all t E M>o, 
it holds that 



( 24 ) 
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a 

and, with da := (— 1 )® IT (OJ + 0 ^ introduced in 

i=Q 


U^^\t) 


nf 

Mi + IW0f)“ 



“ (i + |W0f)“+“ 


• ( 25 ) 


Proof of Assertion\3\ Expression (l2^ is readily obtained from the general Leib¬ 
niz rule. In order to establish (1251) . let / : M>o —)■ M>o be defined as f{z) := 
The feedbaek law ean then be rewritten as U{t) = —j5b^x{t)[foG]{t). 
Using the general Leibniz rule we get that, for any k eN>i and any t G M>o, 
U^k\t) = -/ 3 Z 7 V*)(t)[/oG](t)-fEf^i {^^[foGp^x^^-^Xt). Thanks to Laa 
Di Bruno’s formula (Lemma [T]), we obtain that, for eaeh I G |l,k] and eaeh t G 
M>o, [/oG]W(t)=l'^i/“)(G(t))B;,«(G(i)(t),...,G('-«+i)(0).Sinee/(«)(z) = 
(~l)“n“=o(of+ for all a G |l,k]. Assertion |3] is proven. □ 


We now proceed with Step 2. Set R := min [Rq, ... . We now prove by 

induction on j G |0,p] that there exist /3y > 0 such that, for any /3 < /3j, each 
trajectory of x =Ax + bv{x) satisfies sup{|G*^“'i^(t)| : t > 0} < ^ for ji G |1,7]. 

Note that this in turn ensures that the feedback law dH) is a y-bounded feedback 

law by {Rj,)o<j,<j. 

We start by j = 0. Since sup{|v(x)| : x G M”} < /3, the base case follows by 
choosing /3o = R. Now assume that, for a given j E |0,/? — 1], there exists Pj > 0 
such that for any /3 < /3j the feedback law dH) is a y-bounded feedback law by 
{Rji)o<ji<j- Using Assertion^ we get that for each /3 < /3/ there exists for each 
i G |l,n] and k E |1,7+1] a multivariate polynomial Pi^k : ^ M of degree 1 

(which not depend on /3) such that 


kf^(0l l•*^2(0l): vt>0. (26) 

Lrom (l24l) and (l26l) it follows that, for each k E |l, 7 -t-l], there exists a multivari¬ 
ate polynomial : M" —)■ M of degree 2 , which not depend on /3, such that 


gW (t) I < L,( |xi (01 ,..., 1^2(01 ), Vt > 0. (27) 


In view of Remark 01 (l26l) . and (1271) . we conclude that, for each / G |1,7-|- 1] 
and each a E |1, /], there exists a multivariate polynomial Pi a : M” —)■ M of degree 
2 a -|- 1 , which not depend on /3, such that, for any t > 0 , 

|B;,„(G(i)( 0 ,...,G('-«+i)( 0 )(>^x(^-')( 0 l<n«(|u( 0 l,---,k( 0 l)- (28) 
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Since Pi^a and Pij^i are respeetively of degree 2a + 1 and 1 and reealling that 
a> 1/2, we eonelude that, for eaeh I e |1, y + 1] and eaeh a E |1, /], there exists 
a positive eonstant M; ^ ^ sueh that 

sup{|P,,,(v)/(l + ll^f )“+1 : X e M"} < (29) 

and, for eaeh i E |l,n], there exists a positive eonstant Qij sueh that 

sup{|P,j+i(x)/(l + ||xf)«| : X G M"} < Qi,j. (30) 

Let /3 < 15j and x{-) be a trajeetory of x = Ax+ bv{x). Thanks to Assertion [3] 
(with k = j + 1), we get that, for all t E M>o, \ with 




Ml + IW')lh“ 1-1 V I )h° (i + |W<)lh“+“ 


where Bi^aiG{t)) = Bi^a{G^^\t),.. Therefore a straightforward 

eomputation using (|2^ . (ITTI) . (|2^ . (|29|) . and (l30l) leads to the existenee of a pos¬ 
itive eonstant Mj sueh that |t/*''^+^)(t) | < [5Mj, for all positive time. Thus for all 
[5 < mm{[5j,R/Mj}, it follows that sup{|t/(-^+^)(t)| : t > 0} < ^. This ends the 
induetion on j and eoneludes the proof of Theorem [2l 


4. Numerical examples 

4.1. The triple integrator 

In this subseetion, we illustrate the applieability and the performanee of the 
feedbaek law proposed for Case 1 on a partieular example. We use the proce¬ 
dure deseribed in Seetion 13.11 in order to eompute a 2-bounded feedbaek law by 
(2,20,18) for the multiple integrator of length three. Our set of saturation fune- 
tions is C7i = a 2 = <73 = <7, where o is an saturation funetion with eon- 

stants (2,1,2,1) given by r if |r| < 1, sign(r)(—4-1- 15|r| — 18r^ - 1 - 10|rp — 2r‘^) 
when 1 < |r| < 1.5, 2sign(r)(25 — 60|r| -f- 54r^ — 21|rp -|- 3r^) if 1.5 < |r| < 2, 
and 2sign(r) otherwise. We ehoose = 2/5, = 1/5, = 1/12, and 

= 1 /24. Following the proeedure, we obtain that the two first time derivatives 
of the eontrol signal t/(-) =T(y(-)) = T(//x(-)), with// = (/?i,h 2 ,/? 3 ) where hi = 
( 1 /A^, 2 /A, 1 )^, ^2 = ( 0 , 1 /A, 1 )^ and ^3 = ^ 3 , satisfy sup{|M*^^^(t)| : t > 0 } < 
7.91/A2 + 4.35/A and sup{|m( 2 )(^)| : t > 0} < 26.2/A + 396/A2+ 1147.2/A3 + 
125.2/A4. 
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One get that sup{|t/(^)(t)| : t > 0} < 0.9, and sup{|[/(^)(t)| : t > 0} < 18 by 
ehoosing A =6.5. Observing that the amplitude of U is below 2 by eonstruetion, 
this eonfirms the faet that this is a 2-bounded feedbaek by (2,20,18). The desired 
feedbaek is then given by 

V (;c) = - C7 (x3 + ^ a (5 + ^C7 (24(h[x))))) ). 

This feedbaek law is tested in simulations and the results are presented in Fig¬ 
ure [2l Trajeetories of triple integrator with the above feedbaek are plotted in grey 
for several initial eonditions. The eorresponding values of the eontrol law and its 
time derivatives up to order 2 are shown in Figure [2l These grey eurves validate 
the faet that asymptotie stability is reached and that the control feedback magni¬ 
tude, and two first derivatives, never overpass the prescribed values (2,20,18). In 
order to illustrate the behaviour of one particular trajectory, the specific simulation 
obtained for initial condition xiq = 446.7937, X 2 q = —69.875 and a: 3 o = 11.05 is 
highlighted in bold black. It can be seen from Figure [2] that our procedure shows 
some conservativeness as the amplitude of the second derivative of the feedback 
never exceeds the value 2 , although maximum value of 18 was tolerated. 
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Figure 2: Evolution of the states, the control and its derivative up to order 2 for a set of initial 
conditions. 
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4.2. The harmonic oscillator 

We finally test the performance of the control law proposed for Case 2 through 
example of a 1-bounded feedback law by (2,2) for an harmonic oscillator. We 
consider the following system ii = 5x2, -^2 = —5xi -1- u. 

In accordance with Theorem^ we take u{x\^X 2 ) = —/3x2 /\J 1 +xj+X 2 with 

P = (—5 -f \/4T)/4. The behaviour of the resulting closed-loop system and the 
corresponding values of the feedback and its first time derivative are shown in 
Figure[3]for initial conditions xio = 2 and X 2 o = —2. It can be seen that the values 
of the control and its time derivative stay below 2 as desired. 




Figure 3: Evolution of the states, the control and its first derivative for initial conditions xio = 2 
andx2o = —2. 
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